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Abstract 

We perform detrending moving average analysis (DMA) and detrended fluctuation analysis (DFA) of the WTI crude oil futures 
prices (1983-2012) to investigate its efficiency. We further put forward a strict statistical test in the spirit of bootstrapping to verify 
the weak-form market efficiency hypothesis by employing the DMA (or DFA) exponent as the statistic. We verify the weak-form 
efficiency of the crude oil futures market when the whole period is considered. When we break the whole series into three sub-series 
separated by the outbreaks of the Gulf War and the Iraq War, our statistical tests uncover that only the Gulf War has the impact of 
reducing the efficiency of the crude oil market. If we split the whole time series into two sub-series based on the signing date of the 
North American Free Trade Agreement, we find that the market is inefficient in the sub-periods during which the Gulf War broke 
out. We also perform the same analysis on short time series in moving windows and find that the market is inefficient only when 
some turbulent events occur, such as the oil price crash in 1985, the Gulf war, and the oil price crash in 2008. Our analysis may 
offer a new understanding of the efficiency of the crude oil futures market and shed new lights on the investigation of the efficiency 
in other financial markets. 

Keywords: Crude oil futures, Weak-form efficiency, Detrended fluctuation analysis, Detrending moving average, Bootstrapping, 
Econophysics 
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1. Introduction 

The unfolding financial crisis triggered by the U.S. sub- 
prime mortgage crisis in 2007 has been viewed as the most 
severe in the past centenary after the U.S. stock market crash 
in 1929. The latest financial crisis is followed by the global 
economic recession and the European sovereign debt crisis. 
During this period, we have witnessed the boom and bust 



of many financial and comm odities bubbles dSornettd 12009 



^ ISornette and Woordardi l2010t). The 2006-2008 crude oil bub- 



ble was a significant example (ISornette et all 120091) . The West 



Texas Intermediate (WTI) oil prices experienced an accelerated 
rise followed by a spectacular crash in July 2008. When mea- 
sured in terms of daily close prices, the peak of USD 145.29 
per barrel was reached on July 3, 2008 and a significant low of 
USD 33.87 was seen on December 19, 2008, which is a level 
not seen since 2004. The historical high was USD 147.2 per 
barrel, recorded in the middle of the day of July 11, 200 8. 



In May 2008 before the July peak, ISornette et ail J2009h 



utilized the l og-period i c powe r-law (LPPL) models for finan- 
cial bubbles ( Sornette , 2003alfb ) to analyze the oil prices and 
confirmed the presence of a speculative bubble. Furthermore, 
they successfully predicted the crash of the oil bubble with 
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high precision. The original prediction was posted on the 
arXiv (0806 1170v l) on 6 June 2008. A few weeks later, 



Drozdz et alJ 42008) successfully predicted the Brent oil price 



crash (see ar Xiv:0802.4043y2 , 24 June 2008) based on the 
LPPL models. lYu et alJ (120081) proposed a neural network en- 
semble learning paradigm based on empirical mode decompo- 
sition to forecast crude oil prices and successfully predicted the 
turning point and the oil price plummet. More generally, there 
are a lot of studies showing that the oil prices are predictable 
both in turbulent periods and in normal phases. For instance, 



Sha mbora and Ross iter ( 2007) used an artificial neural network 



model with moving average crossover inputs to predict prices 
in the crude oil futures marke t and u ncovered significant prof- 
itability, and Wang and Yand ( 2010h found significant out-of- 
sample predictability of intraday high-frequency (30-min) WTI 
crude oil futures prices in a bear phase (December 1 , 2000 to 
November 31, 2001) and a bull phase (November 1, 2006 to 
October 3 1,2007). 

A direct consequence of the predictability of oil prices is that 
the oil futures market is inefficient in the weak form. In addi- 
tion, the formation of bubbles and its bust imply that the degree 
of market inefficiency may evolve in a dynamic way. Indeed, 
quite a few studies have been conducted to directly test the ef- 
ficiency of the oil markets, most of which devoted to estimate 
the Hurst index H using different techniques based on the fact 
that a process with the Hurst in dex significantly differing 0.5 



does not execute a random walk (Mandelb rot and Wallislll969 



Mandelbrotl 1983 , 1997 ). These investigations of oil market 
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efficiency can be classified into two related categories. Some 
researchers studied oil prices in relatively long periods, while 
others put emphasis on the evolving behavior of market effi- 
ciency in moving windows of relatively short length^. In the 
rest of this section, we will review the literature along these two 
lines and provide some discussions on existing works. 

1.1. Crude oil market efficiency over long period 



the Hurst index, the larg er the curvature) (iBrvce and Sprague , 



Alvarez-Ramirez et alj (|2002j) performed R/S analysis on 



WTI and Brent crude oil prices over the period from November 
1981 to March 2002. They found that H = 0.567 for WTI crude 
oil prices and H = 0.585 for Brent crude oil prices. Using the 
height-height correlation method (also called structure function 
approach in turbulence), they found that H « 0.47 within the 
sc aling range |T, 181. 



Serletis and Andreadisl (120041) performed R/S analysis on the 



WTI crude oil prices over the period from 2 January 1990 to 
28 February 2001 and found that H = 0.85, which is a sur- 
prisingly large number indicating very strong autocorrelations 
(predictability) in the price time series. When they applied the 
spectral analysis on the prices in the frequency domain, the scal- 
ing exponent was found to be a — 2.03, which corresponds to 
H = (1 + g)/2 = 1.52. It means that the re turn time series 
has H = 0.52. ISerletis and Andreadisl (120041) also calculated 



the structure function of the price time series and found that 

H = 0.69. 

Serletis and Rosenberg! ( 2007 ) performed de trending moving 



average (DMA) analysis on the (logarithmic) WTI crude oil fu- 
tures prices over the period from July 2, 1990 to November 1, 
2006 and reported that H as 0.50. This result indicates nicely 
that, on average, the WTI market was efficient in the weak form 
over the sample period. We also notice that the DMA fluctu- 
ation function has a better power-law scaling than most other 
H urst index estimato rs review ed in this work. 

Elder and Serletis! (12008b adopted a semi-parametric 



wavelet-based estimator to estimate the fractional integration 
parameter d of the WTI oil prices sampled over the period 
from January 3, 1994 to June 30, 2005. They found that the 
fractional integration parameter was d = -0.162 when the 
Daubechies-4 wavelet was adopted and d = -0.243 when the 
Daubechies-6 wavelet was used. Since H = 0.5 + d, these 
results indicate that the WTI oil prices were anti-persistent 
o ver the sample period. 



Alvarez-Ramirez et al.1 (|2008) applied the detrended fluctu- 



ation analysis (DFA) to WTI crude oil prices for the sample 
period 1987-2007. The DFA fluctuation function exhibited a 
crossover at s x = 25 business days separating two power laws 
with very nice scaling behaviors. They found that the DFA 
scaling H = 0.617 + 0.005 for the short time horizons and 
H = 0.507 + 0.004 for the long time horizons. Because the 
DFA fluctuation functions of fractional Brownian motions al- 
ways exhibit downward curvatures at small scales (the smaller 



'it is interesting to note that|Green and Mq5| ll99lh investigated the spot 
and "official" OPEC prices of crude oil for the sample period 1978-1985 and 
found that efficiency was rejected for the sample as a whole but with an im- 
provement over time. 



20121: IShao et all 120121) . we argue that the scaling relation 
with H = 0.507 + 0.004 at large scales reflects the intrinsic 
behavior of the oil prices. It is intri guing to point out that 



this re sult is almost identical to that of ISerletis and Rosenberg 



d2007h. In contrast, with a lo nger sample from 1986 to 2009, 
Alvarez-Ramirez et al. (2010) reported a different DFA fluctua- 



tion function that could not be fitted by two power laws; rather, 
th e function showed multisc aling behaviors. 

Charles and Darne d2009h conducted the variance ratio test 



on the WTI and Brent oil markets over the period 1982-2008. 
They divided the whole samples into two sub-periods at the 
end of 199fl They found that the Brent market were effi- 
cient on both sub-periods, while the WTI market was efficient 
on the 1982-1993 sub-period but inefficient on the 1994-2008 

su b-period. 

Gu et al. ( 2010h also adopted the detrended fluctuation anal- 



ysis to the WTI and Brent oil prices sampled from 20 May 
1987 to 30 September 2008. They also observed a crossover 
in each DFA fluctuation function at s x = 28 for both markets. 
They found that H = 0.53 when s < s x and H = 0.45 when 
s > s x for WTI and H = 0.54 when s < s x and H = 0.47 
when s > s x . We argue that these results hardly support 
th e inefficiency conjecture and roughly c onsistent with th ose 
of I Alvarez-Ramirez et all. d2008l) . Further, iGu et alj d2010h di- 
vided each sample into three time periods separated by the mid- 
dle of the Gulf War (February 24, 1991) and the Iraq War 
(March 30, 2003) and performed DFA on each time period. 
Crossover phenomena were also observed with s x = 21 for the 
first period and s x = 20 for the second and third periods. The 
Hurst indexes were found to vary in [0.52, 0.64] when s < s x 
a nd in [0.36, 0.491 whe n s > s x 



Cuna do et al. l l2010l) utilized different methods to estimate 



the fractional integration parameter d of WTI crude oil prices 
over the period from April 4, 1983 to Se ptembe r 9, 2008, in- 
cluding the modified R/S analysis of Lo ( 19911) . the rescaled 
variance (V/S) analysis of Giraitis et alTi 20031) . and the Whittle 
semiparametric method of Robinson! dl995 ). Their results that 
the hypothesis d = could not be rejected at the 95% confi- 
dence level indicated that there was little or no evidence of long 
m emory in crude oil prices. 



Wang et al.1 ( 1201 11) analyzed the WTI crude oil prices from 



January 2, 1990 to March 9, 2010 using the detrended fluctua- 
tion analysis and identified three scaling regimes. They found 
that H = 0.546 when s < 23, H = 0.472 when 23 < s < 250 
and H = 0.516 when s > 250. We argue that their re sults 



are consistent w ith those of lAlvarez-Ramirez et al.l d2008l) and 
IGu etaljfeoioh . 



iHe and Oianl d2012l) studied the spot prices of the WTI and 
Brent crude oil using the R/ S analysis, the modified R/ S anal- 
ysis and the V/S analysis of ICaiueiro and Tabakl d2005l) based 



2 The separation of the sample at the end of 1993 is rational: In 1993, 
in order to increase the efficiency of the North American energy industry, 
the United States, Canada and Mexico signed the Nort h American Free Trade 
Agreement underpinning the p r ocess of deregulation I Serletis and Andreadis, 
l2004l ; ISerletis and Rangel-Rui j . l2004l) 
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on the work of iGiraitis et al.l (120031) . They estimated that 
H = 0.5406 with the R/S analysis, H = 0.5166 with the modi- 
fied R/S analysis and H = 0.473 1 with the V/S analysis, which 
are respectively very close to the corresponding Hurst indexes 
0.5308, 0.5122 and 0.4806 of the shuffled time series. It means 
that there is no long memory. 

1.2. Evolution of oil market efficiency in moving windows 

Following the idea of estimating the H urst indexes of stock 
marke t s in moving window s of fou r years dCaiueiro and Tabak , 



not surprising that a same method may result in different esti- 
mates of the Hurst index, because different sample periods and 
different scaling ranges have been chosen. 

There are more than ten techniques that have b een invented 



to de t ect long-range co r relatio n s in time se r ies (|Taaqu et al 



19951: belignieres etail l2006t iKantelhardtl l2009t) . such 



as 



the rescaled range (R/S) analysis dHurstL 1 195 lh7 the wavelet 



transform module maxima approach ( IHolschneiden U988 



Muzv et all 1 199 it iBacrv et al 



sis (FA ) dPeng et a 



(DFA) (Peng et al 



20041) . iTabak and Caiueirol (120071) performed the classic R/S y S j s (DMA) ( Alessio et al 



1993), the fluctuation analy- 
, 1992), the detrended fluctuation analysis 
1994. t he det rending moving average anal- 



analysis on the daily Brent and WTI crude oil futures prices 
over the sample period 1983-2004 and found that both markets 
became more efficient and the WTI market was more efficient 
than the Brent market. 



Alvarez-Ramirez et al.l (|20Q8) adopted the detrended fluctu- 



ation analysis to study the WTI and Brent oil prices in moving 
windows. When the size of the moving window is chosen to be 
175 business days, no evident trend was found in the evolving 
Hurst index, which oscillated around H = 0.617. When longer 
moving windows with size of 300 business days were adopted, 
they found that the Hurst indexes decreased from 0.617 to 0.50 
a nd the markets became effi cient . 

Charles and Darne d2009l) investigated the market efficiency 



of WTI and Brent oil (1982-2008) using the variance ratio tests. 
The time period was divided into two parts. Surprisingly, the 
WTI market turned from weak-form efficiency to inefficiency, 
indicating that the deregulation of NATFA did not improve the 
e fficiency on the WTI cr ude oil market as desired. 

Wang and Liu (2010) tested the dynamic efficiency of the 



WTI crude oil market over sample period from July 13, 1990 
to March 6, 2009. The local Hurst indexes were estimated in 
moving windows using multiscale detrended fluctuation anal- 
ysis. They found that the short-term, medium-term and long- 
term behaviors were generally turning into more efficient over 
time. 

Equipped with t he R/S, V/S and modified R/S analyses, 
He and Oianl (12.0121) studied the dynamic efficiency of the WTI 
and Brent crude oil spot markets in moving windows. They 
found that the three Hurst index functions Hrs, Hmrs and 
Hys have similar U-shaped evolution patters and Hrs > 
Hmrs > Hys at ea ch time. Based on numerical experiments, 



He and Oianl d2012l) argued that the R/S analysis overestimates 
the Hurst index, the V/S analysis underestimates the Hurst in- 
dex, and the modified R/S analysis provides better estimates. In 
this sense, we argue that both the Brent and WTI spot markets 
seem to be efficient since the Hurst index evolved within the 
interval [0.45, 0.5 51, which is ex pected to be not significantly 
different from 0.5 dWeronl 120021) . 



1.3. Discussions 

Our brief review has unfolded the evolution of the topic as 
well as the conclusions. Many different methods have been 
adopted, partially owning to the fact that our understanding of 
the performances of different estimators also evolved. It is also 



2002), and their variants. There is 
no consensus which estimators perform best. However, DFA 
and DMA are argued to be "The Met hods of Choice" in deter- 



2012). 



mining the Hurst index of time series dShao et al. 

One possible way to solve the problem is to perform sta- 
tistical tests. For a given estimator, the Hurst index H of 
the original time series is determined. Then, one can shuf- 
fle the time series many times and estimate the Hurst index 
for each shuffled time series. The Hurst index H is compared 
with the average Hurst index (H) of the shuffled time series. 
If the null hypothesis H = (H) cannot be rejected, the orig- 
inal time series can be argued to possess no long-term cor- 
relations. Although there are several well conducted studies 
( WeronL 2002 : Couillard and Davisonl 2005), this approach has 
not been adopted in the majority of previous works. 



2. Data description 

Our analysis is based on the daily prices of the WTI crude 
oil futures, which are freely available at the web site of the U.S. 
Energy Information Administration ( |http : //www . eia . gov/) . 
Our data cover the period from April 4, 1983 to October 2, 
2012. There are totally 7401 data points. FigureQla) plots the 
evolution of the daily prices of WTI crude oil futures. The out- 
breaks of two wars, the Gulf War broken out on August 2, 1990 
and the Iraq War broken out on March 20, 2003, are shown as 
two dashed vertical lines. The solid vertical line corresponds 
to the date January 1, 1994, on which the North American 
Free Trade Agreement are signed by three countries, the United 
States, Canada, and Mexico. 
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Figure 1: (Color online.) (a) Evolution of the daily prices of WTI crude oil 
futures, (b) Evolution of the returns of WTI crude oil futures. The vertical lines 
mark the outbreak of the Gulf War on August 2, 1990, the signing date of the 
NAFTA on January 1, 1994, and the outbreak of the Iraq War on March 20, 
2003. 
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The daily returns are defined as 
r(t) = \nP(t)-\nP(t- 1) 



(1) 



The evolution of the returns of WTI crude oil futures is illus- 
trated in Fig. Hlb)- One can observe that there are many large 
returns after the outbreak of the Gulf war. In sharp contrast 
to this, the fluctuating behavior of the returns almost keeps the 
same as that of the returns before the outbreak of the Iraq war. 
After signing the North American Free Trade Agreement, one 
can observe that the fluctuations of the returns became gradu- 
ally smaller in the following two years. 

3. Methods 

The detren d ed flu ctuation analysis (DFA) was invented 
by IPeng et al.l (119941) originally to investigate the long-range 
dependence in coding and noncoding DNA nucleotide se- 
quences. The pr o perties of DFA have been extensively s tudied 



jHu 



et a 



Xu et al 



2001; Chenetal 



2002, 2005J; IShang et all 120091 



20091: Ma et all l2010t> . The detrending moving 



average (DMA) analysis is based on the m oving average or 
mobile avera ge technique dCarbonel [2 009), which was first 
proposed by IVandewalle and Ausloosl (119981) to estimate the 
Hur st exponent o f self-a ffinity signals and further developed 



by lAlessio et al.l (120021) by considering the second-order 
difference between the original signal and its moving average 
function. The DMA method has been widely applied to the 
analy sis of real-worl d time s eries (ICarbone and Caste lli . 

IVarotsos et all 12005 



20031: ICarbone et al 



ita et al. 

synthetic signals dCarbone and StanlevL I2004t Ku et al 
Serletis , 2008 ; Gu and Zhoui 2010 ). Extensive numerical ex 



Serletis and Rosenberg, 2007; Arianos and Carbone, 2007 



2004bli IVarotsos et al 



Matsushita et all 120071: ISerletis and Rosenberg! 12009b and 



2005 



periments unveil that the performance of the DMA method are 
comparable to the DFA m ethod with slightly different pr i orities 



under different si tuations ( IXu et all 120051: iBashan et all 12008 



Shaoetalll2012l) 



Both methods can be extende d, for instance, to inves 
tigate multifractal time se ries (ICastro e Silva and Mo reira . 



19971: I Weber and Talknerl l200lt iKantelhardt et all 12002 



GuandZhoul |2010[). high-dime n sional f ractal an d mul 



Gu and Zhoul l2006t ICarbonel 12007 



tifractal objects 

iTiirket all 1201 d). and l ong-te r m pow e r-law c ross-correlations 
dPodobnik and StanlevL 120081: IZhoul 120081 IPodobnik etaT 



Jiang and Zhoul 201 lh 



The DFA and DMA share the same framework as described 
below. In the first step, the cumulative summation series y,- is 
determined as follows 



<r», 



1,2, 



,N, 



(2) 



where (r) is the sample mean of the return series. 

The main difference between the DFA and DMA algorithms 
is the determination of the "local trend function" yj, which is 
dependent of the box size s. In the DFA algorith m, the local 
trend y"; is determined by polynomial fits in boxes dPeng et al 



19941) . In the DMA approach, one calc ulates the moving av- 
erage function y* in a moving window (I Arianos and Carbone , 
2007), 



r(j _i X i-e)T 



w»-7 Z 



yi-k, 



(3) 



*=-LCf-i)ej 



where s is the window size, \_x\ is the largest integer not greater 
than x, \x~\ is the smallest integer not smaller than x, and 9 is the 
position parameter with the value varying in the range [0,1]. 
Hence, the moving average function considers \(s — 1)(1 - 9)~\ 
data points in the past and [(s - 1)9] points in the future. We 
consider three special cases in this paper. The first case 6 
refers to the backward moving average dXu et all 120051) . in 
which the moving average function yj is calculated over all the 
past s - 1 data points of the signal. The sec ond case 9 = 0.5 



corresponds to the centered moving average dXu et all 120051) . 
where y, contains half past and half future information in each 
window. The third case 9 — 1 is called the forward moving 
average, where yi considers the trend of s - 1 data points in the 
future. 

The residual sequence e, is obtained by removing the local 
trend function y* from y,: 



e;0) = yt - yi(s). 



(4) 



We can then calculate the overall fluctuation function F(s) as 
follows 



1 N 

[F(s)f = - J] Ms)] 2 

1=1 



(5) 



As the box size s varies in proper scaling ranges, one can deter- 
mine the power law relationship between the overall fluctuation 
function F(s) and the box size s, 



F(s) ~ s H , 



(6) 



where H signifies the DFA or DMA scaling exponent. 

It is necessary to emphasize that the scaling exponent H ob- 
tained from Eq. (0 should be called "DFA exponent" or "DMA 
exponent" rather than "Hurst index" in the rigorous sense. For 
fractional Brownian motions (fBms), the DFA/DMA exponent 
is identical to the Hurst index, which is related to the power 
spectrum exponent rj by r\ — 2H - 1 and to the autocorre- 
lation exponent y by y — 2 — 2H dTalkner and Webei , l200(A 
Heneghan and McDarbvl 20001 Arianos and Carbone , 20071) . 
In addition, ISerinaldil (120101) provided a nice clarification on 
this issue. If the Hurst index of a fractional Gaussian noise is 
H, then the Hurst index of its associated fractional Brownian 
motion is also H. If one replaces P,- for r, in the cumulative 
summation procedure ©, the resultant scaling exponent can- 
not be termed Hurst index. In this case, many estimators give a 
scaling exponent larger than 1 . 

4. Results 

In this section, we provide the analysis results of the mem- 
ory behaviors in the WTI crude oil futures prices by means of 



4 



DMA and DFA approaches. The investigated sample includes 
the whole series, five sub-series, and local samples in moving 
windows. We perform bootstrapping-based statistical tests to 
access if a time series under investigation possesses long-term 
correlations. 

4.1. The whole sample 

To have an overview of the memory behaviors in the WTI 
crude oil futures, we analyze the whole series by means of 
DMA and DFA approaches. Figure [2] illustrates the fluctua- 
tion function F with respect to the box size s as open squares 
for DMA (8 = in panel (a), 8 = 0.5 in panel (b), and 8 = 1 
in panel (c)) and DFA (panel (d)). One can observe very nice 
power-law scaling behaviors between the fluctuation function 
F and the box size s in each panel and the solid lines are the 
best power-law fits to the data points in corresponding scal- 
ing ranges, which gives the exponents H = 0.527 + 0.004 for 
DMA with = 0,// = 0.503 + 0.006 for DMA with 8 = 0.5, 
H = 0.528 +0.004 for DMA with 8 = 1, and H = 0.501 +0.006 
for DFA, respectively. 
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Figure 2: (Color online.) Plots of the fluctuation function F with respect to the 
box size s for the original series and shuffled series, (a) DMA with = 0. (b) 
DMA with 9 = 0.5. (c) DMA with 6 = 1 (d) DFA. 

For comparison, we also perform the same analysis on the 
shuffled series, which is obtained by reshuffling the original 
series to remove possible memory behaviors. The purpose of 
the comparison is to check whether the original series exhibit 
the same memory behaviors as the shuffled series. The conse- 
quence for a randomly chosen realization of the shuffled time 
series is plotted as open circles in each panel as well. One can 
find that the estimated exponents of the shuffled series are com- 
parable to those of the original series, providing an impression 
that the return series of the WTI crude oil futures may exhibit 
same memory features as its shuffled series. 

To distinguish the memory behaviors in original series and 
shuffled series through a strict statistical way, we propose a sta- 
tistical test in the spirit of bootstrapping. The statistic is defines 
as the DMA or DFA exponent H. In order to obtain an ensem- 
ble of the statistic H, which allows us to generate a distribution 



of the statistic, we reshuffle the original series for 10,000 times 
and estimate the DMA or DFA exponent H s for each shuffled 
series. The null hypothesis is that the original series possesses 
the same memory traits as the shuffled series, i.e., H = (H s ), 
where (H s ) is the mean of all H s values. 

Figure [3] illustrates the probability distribution of the es- 
timated exponent H s of the shuffled series for the four ap- 
proaches. The vertical line in each panel corresponds to the 
DMA or DFA exponent of the original series, which locates at 
the central part of the distribution. 




Figure 3: (Color online.) Probability distributions of the DMA and DFA expo- 
nent H s obtained from the shuffled series. The vertical lines correspond to the 
DMA or DFA exponent H estimated from the original series, (a) DMA with 
= 0. (b) DMA with 6 = 0.5. (c) DMA with 8 = 1 (d) DFA. 

The difference in memory behaviors between the original se- 
ries and the shuffled series is quantitatively described by a two- 
tailed p-value, which is defined as 



p = Prob(|//' 5 -<// s >|>|//-<//' 5 >|). 



(7) 



We find that the p-value is 0.998 for DMA with = 0, 0.793 
for DMA with 8 = 0.5, 0.9757 for DMA with 8 = 1, and 0.902 
for DFA, respectively. The results are presented in Table Q] At 
the significance level of 0.01, we cannot reject the null hypoth- 
esis, which offers the strong evidence supporting no memory 
behaviors in return series of the WTI crude oil futures. 

From our analysis we can conclude that the crude oil market 
is efficient when the whole period (1983-2012) is taken into 
consideration. 

4.2. Several sub-samples 

To understand the influences of the important history events, 
which have great shocks on the crude oil market, on the mem- 
ory behaviors in the WIT crude oil futures, we further apply the 
DMA and DFA approaches to five sub-series. Our purpose is 
to find out whether the exogenous events could change the ef- 
ficiency of the oil market in the defined sub-periods. The five 
sub-series are determined as follows. We first divide the whole 
series into three sub-series with two separating dates. One is 
August 2, 1992, on which the Gulf War broke out, and the other 
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Table 1 : Estimated DMA and DFA exponents of the whole series and the five 
sub-series and the results of bootstrapping tests. BDMA, CDMA and FDMA 
are respectively backward, centred and forward DMA. 
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Whole Sub 1 Sub 2 Sub 3 Sub 4 Sub 5 


BDMA 


H 0.527 0.554 0.463 0.590 0.571 0.530 
(H s ) 0.527 0.502 0.502 0.516 0.525 0.522 
p 0.998 0.341 0.394 0.102 0.204 0.778 


CDMA 


H 0.503 0.502 0.401 0.440 0.433 0.482 
(H s ) 0.498 0.496 0.497 0.497 0.497 0.498 
p 0.793 0.860 0.001 0.029 0.010 0.467 


FDMA 


H 0.528 0.557 0.452 0.591 0.564 0.528 
(H s ) 0.527 0.502 0.510 0.516 0.520 0.520 
p 0.976 0.318 0.158 0.098 0.257 0.823 


DFA 


H 0.501 0.481 0.419 0.496 0.439 0.501 
(H s ) 0.503 0.502 0.503 0.503 0.508 0.503 
p 0.902 0.494 0.001 0.817 0.002 0.931 



is March 20, 2003, on which the Iraq War broke out. The three 
sub-series are denoted as sub-series 1, sub-series 2, and sub- 
series 3 based on their chronological orders. We also separate 
the whole series into two sub-series delimited at March 1, 1994, 
on which the North American Free Trade Agreement is signed 
by the United States, Canada and Mexico to increase the effi- 
ciency of the North American energy industry. The two sub- 
series are denoted as sub-series 4 and sub-series 5 according to 
their time sequence. 

Each sub-series is injected into the same DMA and DFA 
analysis procedure as we have done for the whole series. Fig- 
ure 3] plots the fluctuation function F as a function of the box 
size s as open markers for the five sub-series. Again, one can 
see that there are very nice power-law behaviors, which covers 
more than one order of magnitude, for all the five sub-series in 
each panel. The DMA and DFA exponents are estimated by the 
linear fits to In F with respect to In s in the scaling ranges, illus- 
trated as solid lines through the open markers. The estimated 
exponents of the five sub-series are listed in Table Q] 

The proposed statistical test is also performed to verify 
whether the original sub-series and the shuffled sub-series share 
the same memory behaviors. We repeatedly reshuffle the origi- 
nal sub-series and calculate the DMA and DFA exponent of the 
shuffled sub-series till we accumulate an ensemble of 10,000 
exponents H s for each sub-series and each approach. Then the 
corresponding mean value of shuffled exponents (H s ) and p- 
values can be determined, as reported in Table Q] For compari- 
son, we also list the estimated exponent H, the average shuffled 
exponent (H s ), and the p-value of the whole series for each ap- 
proach in the column "Whole" of Table Q] 

At the significance level of 0.01, all the four approaches pro- 
vide the consistent consequences, that the null hypothesis can- 
not be rejected, for sub-series 1, sub-series 3, and sub-series 5. 
This confirms that the crude oil market is efficient during those 
three sub-periods. We also find that the statistical results of 
DMA with 8 = and 9=1 brightly contradict with the results 
of DMA with 9 = 0.5 and DFA for sub-series 2 and sub-series 



10 



10 



10 



(11) 




□ sub 1 
o sub 2 
sub 3 
° sub 4 
a sub 5 



10 



10 



10 



10 




10 



10 



10 



10' 



10 



10 



10 



10 

s 



10 



10 



10 



10 



10 



10 




o sub 1 
o sub 2 
sub 3 
" sub 4 
a sub 5 



10 



10 



10 



10 




10 



10 



10 



10' 



10 



10 



10 



10 



10' 



10 



Figure 4: (Color online. ) Plots of the fluctuation function F with respect to the 
box size s for the five sub-series. The curves have been translated vertically for 
better visibility, (a) DMA with = 0. (b) DMA with 9 = 0.5. (c) DMA with 
$ = 1 (d) DFA. 



4. The DMA with 9 — and 9-1 approaches offers evi- 
dence in favor of no long-term correlations in sub-series 2 and 
sub-series 4, while the DMA with 9 = 0.5 and the DFA favor 
the presence of long-term correlations. From our intuition, this 
contradiction reveals the presence of certain memory behaviors 
in sub-series 2 and sub-series 4. One more thing to note is that 
both sub-periods include the important event of the Gulf war, 
which implies that the Gulf War had reduced the effectiveness 
of the crude oil market to certain degree. 

4.3. Moving windows 

To understand the memory behaviors in local samples, we 
further perform the DMA and DFA analyses on the local series 
within moving windows. The size of moving windows is de- 
fined as 500 and 1000. In each moving window, the original 
series and the 1000 shuffled series are analyzed by the DMA 
and DFA algorithm to compute the exponents H and H s , which 
allows us to proceed our statistical tests. 

When trying to estimated the DMA and DFA exponents 
through linear regressions of In F against In s in moving win- 
dows, we find that the most tough problem is how to determine 
the scaling range in each plot of the fluctuation function F ver- 
sus the box size s. As we know, the best way to locate the 
scaling range is using our naked eyes. However, there are thou- 
sands of figures, in which the scaling range need to be deter- 
mined. We thereby propose an objective strategy to search the 
scaling range between F and s for each moving window. We 
require that each scaling range for fitting must contain at least 
15 data points, which is close to one orde r of magnitude, the 
minimum requirem ent of the scaling range dMalcai et all 1 19971 



Avnir et all I1998D . For convenience, we fix the fitting window 
with a size of 1 5 data points and slide the fitting window from 
the first data point till the last data point. The data points in each 
fitting window are regressed by the linear least-squares estima- 
tion and the fitting window with the smallest fitting residual is 
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regarded as the scaling range, whose slope is recorded as the 
DMA or DFA exponent. For the shuffled series in each moving 
window, we used the same scaling range as their original series 
so that they are comparable to each other. 

Figure[5]illustrates the evolution of the DMA and DFA expo- 
nents in moving windows with a size of 500 (upper panels) and 
1000 (lower panels) for the four approaches. The date in the ab- 
scissa represents the last day of the observations in each moving 
window. The shadow area corresponds to the 2.5% and 97.5% 
quantiles of the DMA and DFA exponents H s estimated from 
the 1000 shuffled series. Although the curve of the exponents 
H exhibits a strongly fluctuating behavior in each penal, most 
of the fluctuations are within the shadow area, which indicates 
that the crude oil market is efficient in most time periods. 

We observe that the evolutionary trajectories for BDMA, 
FDMA, and DFA share remarkable similarity. We further ob- 
serve that there are some spikes, which exceed the 97.5% quan- 
tile, in plots (a), (c), (d), (e), (g), and (h). In our opinion, these 
spikes have connections with some turbulent periods, which are 
related to the oil price crash in 1985, the Gulf War, and the oil 
price crash in 2008. The phenomena further demonstrate that 
the effectiveness of the crude oil market in local periods is re- 
duced when the market is in a turbulent state. In these plots, 
we also find very low H values. Most of them are within the 
2.5/97.5% quantile intervals with a few exceptions in which the 
scaling ranges are problematic. 

We also notice that the results of DMA with 6 = 0.5 ex- 
hibit very different behaviors in the period from 1990 to 1993 
in plots (b) and (f), when compared with the results of the other 
three approaches. To find out why this phenomenon occurs, we 
illustrates in Fig. [6] the fluctuation function F(s) as a function 
of the box size s for three typical windows with their ending 
dates on February 1, 1991, February 3, 1992 and February 1, 
1993. Obviously, we can find that there are crossover behaviors 
for all the three curves in both panels. Our strategic can only 
locate the scaling range, shown as solid lines, in the range of 
s < 50, which leads to the underestimation of the DMA expo- 
nent. However, the crossover phenomena are not observed for 
the other three approaches in the same window. If we aban- 
don the proposed criterion for determining the scaling ranges 
and fit the data with s > 50, the estimates of Hurst indexes 
become normal and comparable to other approaches. We note 
that the shape downward spikes in the plots for BDMA, FDMA 
and DFA have the same behaviors. In conclusion, we argue that 
these very small values of H do not provide conclusive evidence 
for the presence of local anti-persistence. 
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Figure 6: (Color online.) Plots of the fluctuation function F with respect to 
the box size s for three windows with their ending date on February 1, 1991, 
February 3, 1992 and February 1, 1993. (a) Window size 500. (b) Window size 
1000. 



tigated period 1983-2012. Our statistical tests have also been 
carried out on five sub-series in two ways. The two sub-series 
covering the Gulf War were found to exhibit long-memory be- 
haviors, which means that the Gulf War had an impact of reduc- 
ing the efficiency of crude oil market. And we also discovered 
that there were no long-term correlations in other three sub- 
series. 

By adopting the metho d of moving windows inspired by 
Caiueiro and Tabakl (120041) . we further conducted statistical 
tests on the observations in each window. This offers us a new 
insight on the understanding of the evolution of DMA and DFA 
exponents. We found that there are a lot of fluctuations in the 
exponent curve; however, most of the fluctuations are inside the 
2.5%-97.5% quantile intervals of the exponents of shuffled se- 
ries. This implies that the crude oil market is efficient in most 
time periods. The excess of the exponents over the 97.5% quan- 
tile of the shuffled series for the DMA and DFA curves can be 
connected to some turbulent events, including the oil price crash 
in 1985, the Gulf war, and the oil price crash in 2008. Our re- 
sults further revealed that the crude oil market requires a certain 
time to digest the information. If the investigation time scale is 
longer than the digesting time, the market looks efficient, and 
if the investigation time scale is shorter than the digesting time, 
the market looks inefficient. 
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5. Conclusion 

In this paper, we performed detailed analyses on the weak- 
form efficiency of the WTI crude oil futures market by means 
of DMA and DFA approaches. By using the DMA and DFA 
exponents as the statistic, we proposed a strict statistical test in 
the spirit of bootstrapping to check whether the original series 
exhibit the same memory behaviors as the shuffled series. 

By performing the statistical tests on the whole series, we 
demonstrated the efficiency of the crude oil market in the inves- 
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